The Connection Between Density Matrix Method, Supersymmetric Quantum 
Mechanics and Lewis-Riesenfeld Invariant Theory* 
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This paper is concerned with the connection between density matrix method, supersymmetric 
quantum mechanics and Lewis-Riesenfeld invariant theory. It is shown that these three formula- 
tions share the common mathematical structure: specifically, all of them have the invariant oper- 
ators which satisfies the Liouville-Von Neumann equation and the solutions to the time-dependent 
Schrodinger equation and/or Schrodinger eigenvalue equation can be constructed in terms of the 
eigenstates of the invariants. 
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I. LEWIS-RIESENFELD INVARIANT THEORY 

Lewis-Riesenfeld invariant theory can be applied to the solutions of time-dependent Schrodinger equation. 

In order to illustrate the Lewis-Riesenfeld invariant theory [1] easily, we consider a one-dimensional system whose 
Hamiltonian H(t) is time-dependent. According to the Lewis-Riesenfeld invariant theory [1], a Hermitian operator 
I(t) is called invariant if it satisfies the following invariant equation, i.e., the Liouville-Von Neumann equation (in the 
unit h ^ \) 



The eigenvalue equation of the time-dependent invariant |A„,t) is given 

/(t)|A„,t) = A„|A„,t), 

where 

0. 



dXn 
dt 



The time-dependent Schrodinger equation (in the unit fi = 1) for the system is 



dt 



H{t)\^{t))^. 



(LI) 



(1-2) 



(1-3) 



(1.4) 



In terms of the Lewis-Riesenfeld invariant theory, the particular solution |A„,i)j of Eq.(1.4) differs from the eigen- 
function |A„,t) of the invariant I{t) only by a phase factor exp[i0„(i)], then the general solution of the Schrodinger 
equation (1.4) can be written as 



l*W>. = ECnexp[-0„(i)]|A„,t), 

n. ^ 



(1.5) 



where 



An, t 



H{t 



dt' 



K,t ) dt , 



C„ = {Xn,t = 0\^it^O))^ 



(1.6) 



*This paper was initiated three years ago (in 2000). 
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l^m t) g = exp[i(f>n{t)] |A„, t) (n = 1, 2, • • •) are said to form a complete set of the solutions of Eq.(1.4). The statement 
outlined above is the basic content of the Lewis- Riesenfeld invariant theory. 

A time-dependent unitary transformation operator can be constructed to transform I{t) into a time-independent 
invariant ly = V^t)I{t)V{t) [2,3] with 

Iv = A„ |A„) , (1.7) 

\Xn) =VHt)\Xn,t). (1.8) 

Under the unitary transformation V{t), the Hamiltonian H{t) is correspondingly changed into Hv{t) 

Hvit) = V\t)H{t)V{t) - V\t)z^. (1.9) 

In accordance with this unitary transformation method [2] , it is very easy to veriiy that the particular solution |A„, t)^Q 
of the time-dependent Schrodinger equation associated with Hv{t) 

i^-^^^=Hy{t)\X^,t)g, (1.10) 

is different from the eigenfunction |A„) of ly only by the same phase factor exp[j(/)„(f)] as that in Eq.(1.5), i.e., 

\Xn,t) sQ = ^'^'P\r:(t^n{t)\\Xn) . (1.11) 

Substitution of |A„,f)^Q of Eq.(l.lO) into Eq.(l.ll) yields 

</.(t)|A„) =iIy(t)|A„), (1.12) 

which means that Hv{t) differs from Iv{t) only by a timc-dcpcndcnt multiplying c- number factor. It can be seen from 
Eq.(1.12) that the particular solution of Eq.(l.lO) can be easily obtained by calculating the phase from Eq.(1.12). 
Thus, one is led to the conclusion that if the V{t), ly, Hv{t) and the eigenfunction |A„) of ly have been found, the 
problem of solving the complicated time-dependent Schrodinger equation (1.4) reduces to that of solving the much 
simplified equation (1.10). 

II. DENSITY MATRIX METHOD 

Density matrix method has many helpful applications to the semiclassical theory of laser [4] . 

The density matrix operator is defined to be p{t) = |\l/(f))(\E'(t)|, the wavefunction \^{t)) of which agrees with 
the time-dependent Schrodinger equation. It follows that the density operator satisfies the following time-evolution 
equation 

^ + i[p(f),if(t)]=0, (2.1) 

which is just the Liouville-Von Neumann equation (1.1). This, therefore,means that the density operator p{t) is the 
Lewis- Riesenfeld invariant. 

In what follows we will discuss the density matrix method that solves the time-dependent Schrodinger equation 
governing the semiclassical laser process. As an illustrative example, we consider a two-level atomic system (|a) and 
\b)) interacting with the classical Maxwellian electromagnetic fields. The time-dependent wavefunction \'^{t)) can be 
expressed in terms of \a) and |6), namely, 

|*(i))=c„(t)|a)+C6(t)|6), (2.2) 
and the density matrix can be rewritten as 

/^(i) = il:) « 4) = { '?$ l:i ) (2-3) 

in the representation of \a) and |6). 



Ch ) ^ \ CbCl CbCl 
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The Hamiltonian of two-level atomic system takes the form 

H = Ho + V, 

where the atomic free Hamiltonian Hq and the interaction Hamiltonian V are respectively of the form 



Ho 



OJa 

(Jb 



V 



Vab 
Vba 



Substitution of (2.3)-(2.5) into (2.1) yields 

p'aa = - [iVabpba + C.c], p'bb = iVabpba + C.C, p'ab = -i{uJa " i^b)Pab + iVab{Paa " Pbb) 



(2.4) 



(2.5) 



(2.6) 



with paa = CaC*ai Pab = Cacl, Pba = Cf,c* and pbb = Cbcl ■ So the rcsolutiou of the Schrodinger equation can be 
ascribed to the solution of the Eq.(2.6) {i.e., the Liouvillc-Von Neumann equation (1.1)), which is similar to the case 
in Lewis- Riesenfeld invariant theory, where we should obtain the expression for the invariant I{t) via the Liouville-Von 
Neumann equation and solve the eigenvalue equation (1.2) of I{t). 



III. SUPERSYMMETRIC QUANTUM MECHANICS 



Supersymmetric quantum mechanics is applicable to the energy eigenvalue problems of one-dimensional potential 
walls. 

The sequence of ideas of supersymmetric quantum mechanics is to obtain the eigenstates of the following two 
Hamiltonians 



via the eigenstates oi H = ^ 



+ V{x), where V±{x) agrees with the following two equations: 



\[V+{x) + V.{x)\ = W{xf, \[V+{x) - V.{x)] 



'2m 



W'{x) 



(3.1) 



(3.2) 



with W'{x) denoting the derivative of W{x) with respect to x, and W{x) is so defined that the ground state 'il)o{x) of 
H is 



'2m 



tpo{x) = exp 

Define the lowering and raising operators A and as follows 

A = 



dxW{x) 



h d , .+ ?i d 

+W{x), AT = 

V2mdx 



/2m dx 



W{x), 



(3.3) 



(3.4) 



and therefore wc have = A'' A and = AA^ and the commuting relation [A, A^ 



2h 



W'(x). 



By using the properties of lowering and raising operators A and , we can obtain the eigenstates of if±, which 
are constructed in terms of the eigenstates of H. 

It is readily verified that here the Hamiltonian H serves as the invariant that satisfies the Liouville-Von Neumann 
equation. Since here H is time-independent, we should calculate the commuting relation [H^,H] only. According to 
Eqs.(3.2), we obtain V+ = W{xf + --^W'{x) and consequently 



V+ = 



h 



2 d'^ipajx) 



2m ipo{x) 



(3.5) 



In accordance with the eigenvalue equation [— + V{x)]''Po{x) = So^o(x), we get V+(a;) = V{x) — eq. So 
[H+,H] = [V+(x) — y{x),—-^-^] = 0, because of V+{x) — V{x) = —Sq. This, therefore, means that the Hamiltonian 
H is an invariant. 

For the detailed applications of supersymmetric quantum mechanics to the energy eigenvalue problems of one- 
dimensional potential walls and related problems, Readers may be referred to the references [5 8]. 
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IV. CONCLUDING REMARKS 



In general, the Lewis-Riesenfeld invariant theory is appHcable to the multi-level atomic quantum system interacting 
with the second-quantized electromagnetic fields [9,10], while the density matrix method is appropriate to treat the 
multi- level atomic quantum system interacting with the classical Maxwellian electromagnetic fields (waves) [4]. The 
former formulation can also be applied to the motion of a charged spinning particle in a classical magnetic fieldcite 
[11] or of a photon moving in a sufficiently perfect curved fibercite [9]. 

In the present paper we show that the above three formulations share the common mathematical structure: specif- 
ically, all of them have the invariant operators which satisfies the Liouville-Von Neumann equation and the solutions 
to the time-dependent Schrodinger equation and/or Schrodinger eigenvalue equation can be constructed in terms of 
the eigenstates of the invariants. 
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